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Abstract 

We use the methods of group theory to reduce the equations of motion of the 
CP^ model in (2+1) dimensions to sets of two coupled ordinary differential 
equations. We decouple and solve many of these equations in terms of ele- 
mentary functions, elliptic functions and Painleve transcendents. Some of the 
reduced equations do not have the Painleve property thus indicating that the 
model is not integrable, while it still posesses many properties of integrable 
systems (such as stable "numerical" solitons). 

I. INTRODUCTION 

Over the past few years, it has become clear that many physical processes can be de- 
scribed well in terms of various partial nonlinear differential equations. The areas providing 
such equations, range from solid state physics, hydrodynamics and particle physics to bio- 
physics and biochemistry. As the equations are nonlinear, in general, they are hard to solve; 
in fact, so far no general method of solving these equations is known and each equation has 
to be treated on its own. However, a particular class of equations, which are derived from 
the so-called integrable models can be solved using some very general techniques |l[| (inverse 



scattering methods, Bdcklund transformations, ... ). These equations are very special and 
their solutions have very special properties. Many of these equations have solutions which 
are localised in space and propagate at a constant speed. Such solutions, usually called 
solitons, or extended structures in general, have received a lot of attention in recent years. 
However, most of these equations depend on two space (or space-time) variables and, as 
such, can only describe phenomena which are quasi-twodimensional. When they involve 
more variables, either all variables come in a very nonsymmetric way or the models are very 
special. 

Most applications in nature involve 3 spatial dimensions, and in many applications all 
spatial variables come on an equal footing. When the applications involve, for example, 
particle physics or relativity, the underlying models are Lorentz covariant. Such models 
are, generally, nonintegrable and the methods mentioned above do not apply. On the other 
hand, some of them can be studied numerically. Such studies have involved full simulations 
of similar models in (2+1) dimensions or simulations of various reduced models (i.e. approx- 
imations to the original models). Some of these studies [Q] have found that even though 
the models were not integrable the behaviour of their extended structures resembled the 
behaviour that was expected had the models been integrable (i.e. the structures preserved 
their shapes well, there was little radiation etc.) Morever, the approximate methods |^ gave 
results which were virtually indistinguishable from the results of full simulations. Hence 
we feel that the behaviour of integrable models may not be that very unusual; other mod- 
els which, strictly speaking, are not integrable maybe be "almost" so; with their extended 
structures showing in their behaviour very little difference from what may be expected had 
the models been integrable. Moreover, various approximate methods work well and may be 
used to provide some insight to the behaviour of the solutions of the full equations. 

Most of the results, which involve models in (2+1) dimensions, were obtained in the 
so-called CP^ model (also called the S*^ or 0(3) model) and its modifications. |0 This 
model, in its original version [Q, is probably the simplest model in (2+1) dimensions which 
is relativistically covariant and which admits the existence of localised soliton-like solutions. 
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The CP^ a model is defined by the Lagrangian density 

= f (9^0) ■ (dj), (1.1) 

together with the constraint (p .(p = 1. The Euler-Lagrange equations derived from ( p. .11) are 

d>'dj+id''$-dj)$=6. (1.2) 

The model can be modified by the addition of further terms to the Lagrangian den- 
sity. The terms that have been studied the most extensively involve the (2+1) dimensional 
analogue of the "Skyrme" term and various "potential" terms [0,^. They were added, 
primarily, to stabilise the soliton-like structures. In the original model the soliton-like struc- 
tures were not really stable; any perturbation would induce their shrinking, or expanding 
which they could do without any cost of energy due to the conformal invariance of the 
pure CP^ model. Apart from curing the shrinking, and inducing also weak forces between 
soliton-like structures the additional forces had little effect on the dynamics of these struc- 
tures. Moreover, the affects of nonintegrability of these models, were also not that different 
from similar effects in the pure CP^ case. Hence the dynamics of these models was de- 
scribed well by the dynamics of the CP^ case. The same was true when one looked at the 
approximate methods 

These observations suggest that a lot can be learnt from looking at exact solutions of 
the CP^ model using the group theoretical method of symmetry reduction ||8|-[T0|. This 



method exploits the symmetry of the original equations to find solutions invariant under 
some subgroup of the symmetry group (the classic example one can give here involves seeking 
solutions in two dimensions which are rotationally invariant). The method puts all such 
attempts on a unified footing and it has been apphed with success to many equations. The 
method leads to equations whose solutions represent specific solutions of the full equations; 
the solutions are determined locally and the method does not tell us whether these solutions 
are stable or not with respect to any perturbations. 

In our case we would like to apply this method to looking for solutions of the original 
CP^ model; from the remarks made above we can hope that these solutions will be also 



approximate solutions of the modified models. Their stability is harder to predict; but 
again, guided by the experience from the numerical simulations we hope that, at least, some 
of them will be stable with respect to small perturbations. 

In order to perform the symmetry reductions of the pure CP^ model in (2+1) dimensions 
we have to decide what variable to use. To avoid having to use the constrained variables 
(0) it is convenient to use the W formulation of the model which involved the stereographic 

— * — * 

projection of the sphere 0-0 = 1 onto the complex plane. In this formulation instead of using 
the (f) fields we express all the dependence on in terms of their stereographic projection 
onto the complex plane W . The fields are then related to W by 

In this formulation the Lagrangian density becomes 

^ - {i + \w\^r ^^-^^ 

where * denotes complex conjugation. 

To perform our analysis it is convenient to use the polar version of the W variables; i.e. 
to put W = Rexp{iip) and then study the equations for R and ijj. The advantage of this 
approach is that the equations become simple; the disadvantage comes from having to pay 
attention that R is real and t/j should be periodic with a period of 27r. (If the period is not 
27r then the solution may become multi- valued etc.) Thus if we find solutions that do not 
obey these restrictions, then these solutions, however interesting they may be, cannot in 
general be treated as solutions of the original model. 

The equations for R and t/j take the form 

dtt'tp - d,,7P - dyyiP + 2 ^|i~^^2) (W^ - d,7j;d,R - dy^PdyR) = (1.5) 

and 

R(l — R^) 

dttR - da^xR - dyyR - (^^^2) (^tMV^ - d^ipdxil^ - dy^dyip) 

2R 



j^^^^iidtRf - {d.Rf - {dyRf) = 0. (1.6) 
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Note, that if we put R = 1 the second equation is automatically satisfied and the first 
one reduces to 



dttip - dxxip - dyyijj = (1.7) 

i.e. the linear wave equation for the phase ip- 

In section |l| we determine the symmetry group of equations (|1.5|) and (1.6). In the 
following section we present coupled pairs of reduced ordinary differential equations (ODE's) 
for all two-dimensional subgroups of the symmetry group. Sections and are devoted 



to the presentation of explicit solutions. We finish the paper with a short discussion of the 
derived solutions, their relation to the solutions known before and their physical relevance. 

II. THE SYMMETRY GROUP AND ITS SUBGROUPS 



The symmetry group of the system (|1.5|) and (1.6) can be calculated using standard meth- 



ods p-[TO||. We actually made use of a MACSYMA package [|TT| that provides a simplified 
and partially solved set of determining equations. 

Solving those we find that the symmetry group has the structure of a direct product, 
namely, 

G ~ SIM(2, 1) ® ^f/(2), (2.1) 

where SIM(2, 1) is the similitude group of (2+1) dimensional Minkowski space (the Poincare 
group extended by dilations). The group SU{2) rotates the components of the fields amongst 
each other. 

The corresponding Lie algebras sim(2, 1) and 5^(2) can be represented by vector fields 
acting on R and ip and the space-time coordinates. A suitable basis is given by two Lorentz 
boosts Ki and K2, one rotation L, three translations Pq, Pi and P2, one dilation D and 
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three su{2) generators X, Y and Z. We have 

Ki = -{xdt + tdx), K2 = -{ydt + tdy), L = -xdy + yd^ 
Po = dt, Pi = da-, P2 = dy, D = tdt + xdx + ydy, 



(2.2) 



(2.3) 



z = a 



Our aim is to obtain solutions of eq. (|1.5|) and (1.6 by the method of symmetry reduc- 
tion [PHTOH. In practice we shall require that solutions are invariant under a two-dimensional 



subgroup of the symmetry group G. This will reduce the original partial differential equa- 



tions ( p..5|) , (1.6) to a system of ODEs. Subalgebras of sim(2, 1) were classified in ref. [T^. A 
two-dimensional algebra {A, B} can be either Abelian, [A, B] = 0, or solvable non-Abelian, 
[A, B] = A. Subalgebras of the direct sum sim(2, 1) © 5^(2) can be obtained by the Goursat 
"twist" method. [|T3iP 



The result is the following. There exist 10 parametric classes of Abelian subalgebras, 
represented by 

{A + aZ, B + bZ} (2.4) 
with a,b (z M., Z as in (2.3) and {A, B} equal to one of the following pairs: 

{K,,P2}; {D,K2 + L}; {K^ + L, P^ - Pi}; {P^.^o-^i}; {L,Po}; 
{Pi,P2}; {Po,Pih {D,K,}; {D,Ly, {D - K,,P,- P,}. 

(2.5) 

Further, there exist 15 parametric classes of nonabelian two-dimensional subalgebras repre- 
sented by 

A + cZ,B cgM, (2.6) 
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where {A, B} is one of the pairs: 

{Ku K2 + L}; {D, P2}; {Ki, Pq - ^i}; {D, Pq - ^i}; {Ki + eP^, Pq - Pi}; 
{D + e{K2 + L),Po- Pi}; {D, Pq}; {D + aK,, + L}; {D + ai^i, P2}; 
{D~K, + e(Po - Pi), + L}; + 7^2 + e(Po + P2), Pi}; {D + aL, Pq}; 
{D + aKi, Po - Pi}; {D + Ki + e{Po + Pi), Pq - Pi}; 
{D + ^K,,K2 + L + e{Po + P,)} 

with a G M and e = ±1. 



III. THE REDUCED EQUATIONS 



For each subgroup (|2.4|) and (|2.6|) we find three invariants, ^, P and P, using standard 



methods. PH1Q| In terms of these we express the two functions R and ijj of (|1.5| ) and (1.6) 

as 

P = P(0, ^ = a(x,y,t) + P(0, e = e(a:,2/,t), (3.1) 

where a and ^ are given for each subalgebra in Table | and Derivatives with respect to 
the variable ^ will be denoted by dots. We introduce the two invariant operators A and 
by setting 

A/ = ftt - fxx - fyy, (V/, \/g) = ftQt - fxQx " fyQy (3.2) 

and consider three cases separately. 
1. (VO^ ^ 0. We put 

- V--'-. ^^-K S^/, (3.3) 



(VO^ " g' (VO^ ' (VO^ ' (VO 

with f7 = f7(0, = /i(0> s = s(0 and / = 
For further use we also introduce 

m = h--h-s (3.4) 
9 
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and 



B = 2^—. 3.5 



The two PDE's ( [L.5|) and (1-6) now reduce to 



(3.7) 

For each algebra the functions ( |3.3|) are given in Table |. 

In order to solve the above system we must decouple its two equations. Putting 
F + h = V v/e first rewrite equation (3.6) as 

with m as in eq. (|3.4|) . 

For m = we solve (3.6) and obtain 

F = A^^-^^9{i) -K A = const. (3.9) 



Next we substitute F into equation (3.7) and obtain a second order ODE for R{i) 
1 + g R-^ 1 + R^ 



(3.10) 



For m 7^ eq. (|3.8|) is inhomogeneous. We can still decouple it by putting 



Um 

F=—. h 3.11 

U 

[R^ + lyg 

Using eq. (3.11) we can rewrite (3.7) as a third order ODE for the auxiliary function 
U{^). If we can solve it we obtain R{C,) from (3.12). However, in this paper we restrict 
our attention to the case of m{^) = 0. 
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2- (V^)^ = 0, = 0, (Va, V^) 7^ 0. The reduced equations decouple immediately and 
we have 

(Va,VOF=-i(Vaf. 

For (Va, V^) 7^ we can integrate directly to obtain 

R ^ { 1 r Aa 



(3.14) 



1 r (VaY 

For (Va, VO = we must have also Aa = 0, (VO^ = 0. Then R{^) and F{^) are 
arbitrary functions of ^. In particular this is true for a = 0. 

3- (V^)^ = 0, = h{^) 7^ 0. This case can actually not occur, since these two equations 
are compatible only for h{C,) = 0. 

IV. ANALYSIS OF SECOND ORDER ODE 
A. General Comments 



In order to obtain explicit solutions we need to solve the ODE ( p.lO| ) for the function 



R{^)- This equation is in the class analysed by Painleve and Gambler |]T5|-|T^; namely it is 
of the form 

y = f{y,y,x) (4.i) 

where / is rational in y and y and analytical in x. If this equation has the Painleve property 
(no movable singularities other than poles) then it can be transformed into one of the 
50 standard equations listed e.g. by Ince [l^. The Painleve test [p!8| , p!9| provides us with 



necessary (but not sufficient) conditions for eq. ( [4. 1| ) to have the Painleve property. The 



solution is expanded about an arbitrary point xq in the complex x-plane in a Laurent series 



00 



y(a;) = ^ afcr'^+P, t = x - Xq, (4.2) 

k=0 
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where p is required to be an integer (usually a negative one). The coefficients are obtained 
from a recursion relation of the form 

Pkak = hk{xo,ao,ai,. . .ak-i). (4.3) 

Since ( [4.2|) is supposed to represent a general solution it must depend on two constants; 
Xo and for some nonnegative r, called a resonance value. This occurs if the function Pr 
satifies Pr = 0. Then is arbitrary and we have a consistency condition, the "resonance 
condition" , hr{xo, oq, ai, . . . a^-i) = (which must be satisfied identically in xq). If the above 
conditions are satisfied, the Painleve test is passed and ( [4.2| ) represents a two parameter 
family of formal solutions (locally, within the radius of convergence of the series). 

Turning our attention to eq. ( |3.10| ) we note that the cases A ^ and A = must be 
treated separately. 



B. The case A^O 

The Painleve test applied directly fails immediately since a balance between the most 
singular terms occurs for p = — |, i.e. we have a movable square root branch point. To 
remedy this problem we put 



^(0 = ^J-UiO, 



(4.4) 



and obtain 



U = U' 



1 1 



2U U-1 



u 
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Ujl + U) 

iu-i) ' 



We can now choose a new variable rj to be 



(4.5) 



(4.6) 



and transform eq. ( [4.5| ) into 



1 



2U U-1 



U (U-l) 



(4.7) 
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with B as in ( p.5|) . 

For i?=const this is eq. PXXXVIII hsted e.g. by Ince |T^. It has a first integral K that 
we use to write a first order equation for U : 



AA'U^ + AKU^ + (8A^ -2B- 8K)U' + AKU - AA' 



(4.8) 



Since we have A ^ we can rewrite ( ^.81 ) as 



(Uy = -AA\U - U{){U - U2){U - U,)iU - U,), 



(4.9) 



where the constant roots of the right hand side of ( |4.9| ) Ui, i = 1 . . . 4 satisfy 



UiU^U^Ui = 1 

A"^ 

U1U2 + f/if/s + UiUi + f/af/s + f/2f/4 + U^Ui = -2 



B 2K 

+ 



2^2 A^ 



U1 + U2 + U3 + f/4 



K 



(4.10) 



Eq. ( [4.9| ) has elementary algebraic and trigonometric solutions, as well as solutions which 
resemble solitary waves or kink-like structures (in the symmetry variable 1]) in the case of 



multiple roots. If all Ui are distinct then the solutions of ( [4. 81) involve elliptic functions PO 
Explicit solutions will be presented in the next section. 

If B in eq. ( |4.8|) is not constant we proceed differently. We again introduce a new 
independent variable 



1] = e 



Jgdi 



(4.11) 



and transform eq. ( [4.5| ) into 
1 1 



For 



U = U' 



2U U-1 



1] 1]'^ U g'^T]'^ U — I 



5 = const 



(4.12) 



(4.13) 
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this is the equation for the fifth Painleve transcendent. 

The values m = and B ^ const occur for the cases 1 and 2 from Table | and we actually 
have 

r/ = e, 5 = 2(6^±^), (4.14) 

so b = const requires a = 0. 

Hence we obtain the solutions 

f/(0 = Py(a,/3,7,5;0 

(4.15) 

a = -/5 = 7 = 0, b = 26^ 

for equations describing the group reductions 1 and 2 (of Table |) in the case when a = 0. 

For 6 = in the cases of these two reductions we have B = 2a^ = const and we obtain 
solutions in terms of equation ( ^Tsf ). 



C. The case A = 0, B = const 



The transformation ( [4.4| ) can also be performed for the case A = 0. We use the first 
integral to write our equation as 



U'^ = AK{U -Ui){U -U2) 

B 



Ui,2 = L + l±^LiL + 2), L = —, Kj^O. 
For K = we find a solution immediately; namely we have 



U = f/oe^^^^ B<0. 



(4.16) 



(4.17) 



On the other hand, for A = in eq. (|3.1CI|) , the Painleve expansion (|4.2| ) gives us p = — 1 
for the leading (most singular) power. Hence the transformation ( [4.4|) is not required and 
so we can transform eq. ( [4.2|) directly into one of the standard forms. 

We put 



R=-i 



(4.18) 
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and obtain 



(4.19) 



If B of eq. ( p.5| ) is constant we can choose r] as in eq. ([4 .61), set yU = 1 and obtain the 
equation 



- P BfZ^-l 

z = — + - 



Z 8\ Z 



(4.20) 



This equation can be integrated directly for i? = 0. For i? 7^ it has a first integral K in 
terms of which we obtain a first order ODE 



= -(Z^ - zl){z^ - zl 



(4.21) 



= -K± VK^ - 1. 



We again have elliptic function solutions. They are however not new. Since R is real (and 
nonnegative) we must require that 



Z = e'''^'^\ ioT fi= 1,0 <ri < 271. 



(4.22) 



The relation between the function (7(1]) and U{r]) of eq. (^.1|) is 



sm a 



1 — sin a 

For B T^const we use the variable (^4.101) , again set /x = 1 and reduce eq. ( [4.19|) to 

Z = ---Z + ^^^. 
Z 7] Ag'^rf Z 



(4.23) 



(4.24) 



For 5, defined by ( |4.13|) , being constant (5 = 5^) eq. ( |4.24|) is a special case of the third 



Painleve transcendent Pin(a, 7, /?, 5 : ifj and so we have as a solution of eq. ( 4.24 ) 

5n 5i 



Z = P:n(0,0,|,-|;,). 



(4.25) 



This equation is, however, not new; it can be transformed into solution ( 4.15 ) by making 
use of relations between special cases of Py and Pm. 
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D. Comments on the Painleve analysis and integrability of model 



For A 7^ eq. ( [4 .51 ) passes the Painleve test for any function g{$,) and /i^(^) — /(^). Indeed, 
we find p = —1 in the expansion (|4.2| ). A resonance is obtained for k = 1; the resonance 
condition is satisfied and so the coefficient oi is a free constant (as is xq). 

The Painleve test only checks whether certain necessary conditions are satisfied. If an 
equation of the type ( [4 .11) does actually have the Painleve property (as opposed to merely 
passing the Painleve test), then it can be transformed into a standard form by a Mobius 
transformation (with variable coefficients) 

±1. (4.26) 

Eq. ( |4.5| ) is already, to a large extend, standardized. Indeed, the coefficient of f/^ has 
poles at f/ = 0, 1 and oo. This puts the equation into Ince's class IV and the residues have 
the correct values. Hence we have a = p = 1, (3 = 6 = in eq. ( ^.26| ). The only remaining 
permitted transformation is that of the independent variable. We have shown above that 
eq. (|0| ) can be reduced to the elliptic equation if and only if B in eq. ( ^^) is constant. It 
can be transformed into the equation for the Py transcendent if and only if its coefficients 
satisfy 

^(^)_2(^)<, = 0. (4.27) 

In all other cases the equation ( [4.5|) cannot be transformed into a standard form and hence 
it does not have Painleve property. 

The situation is exactly the same for A = 0. Eq. ( |3.10|) passes the Painleve test and is 
transformed into eq. ( [4.19|) by a Mobius transformation. The coefficient of has poles at 
Z = and Z —>■ oo with the correct residues. Hence only Z{^) a{^)z(ri{^)^ is permitted. 
Eq. ( |4.19| ) is of Ince's type II and has the Painleve properties if and only if eq. ( |4.27|) is 



satisfied. 

Thus we have shown that eq. ( |4.5| ) has the Painleve property if and only if the coefficients 
satisfy B = const, or Eq. ( [4.27|) . 
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V. EXPLICIT SOLUTIONS 



We have reduced the original system (|1.5|) , (1-6) for the function 

W = Re'^ (5.1) 

to one of the pairs of equations { (|3l9|) , (|3.10| )|, or {(3.14), (3.15)}. 

Let us first look at the pair |(pl9|), (|3.1CI| )|. Eq. (3.9) provides -F(0 by a quadrature, 



once eq. (|3.1CI|) is solved. In section ^\ we have further reduced eq. ( p.lOj ). Using eq. (^^ 



we replace equations for R{C,) by equations for U{^), where f/(0 must satisfy U{^) <0. 

As mentioned above, the algebras 1 and 2 of Table | lead to solutions of the form ( [4. 151 ), 
i.e. the fifth Painleve transcendent Py{C,), for a = 0, 6 7^ 0. 

Algebras 1-23 lead to the elliptic function equation (4.9) for m = 0, A 0, B = const 
and to eq. ( |4.1(j| ) form = 0,y4 = 0,i? = const, K 0. In both cases, elementary solutions 



are obtained in the case of multiple roots of the polynomial on the right hand side of the 
equation. Let us discuss in more detail the real nonsingular solutions of these equations. 
The character of the solution depends crucially on the sign of B in the table {B is defined 
in (|3j)). We have 

• 5 > for algebras 

1 (6 = 0, a ^ 0), 10 (a^ + 6V 0), 11 (a ^ 0), 12 (a ^ 0), 
15 (a^ > b^),16 {b ^ 0), 
17 (a = 0, 6 ^ 0, x^ + y^ -t"^ > 0), 18 (a = 0, 6 ^ 0, + - t'^ > 0), 
19 {a = 0,b^O,x'^ + -t^ > 0). 

• B < for algebras 

2 (6 = 0, a ^ 0), 13 (a ^ 0), 14 {ab ^ 0), 
15 (a^ -6^ < 0), 

17 (a = 0, 6 ^ 0, - - y2 > 0), 18 (a = 0, 6 ^ 0, f -x^ -y^ >0), 
19 {a = 0,b^O,t^ - x^ -y^ > 0). 



(5.2) 



(5.3) 
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In all other cases with m = we have B = 

Let us first run through all elementary functions solutions, remembering that the inde- 
pendent variable is rj given in eq. ( [4 .61 ). 

Localized solutions are obtained precisely for the algebras ( [5 .21 ). From eq. ( [4.16| ) (i.e. 
A = 0) we obtain a kink in R{^) where ^ and function /i(^) are read off from Table |. Two 
cases are to be considered. 

L A = 0, L = -2, K <0, B >0: 

The solution is: 

1. 

1 IB f 
R = ±tanh-J-{r]-7]o), F = - J h{r]) cIt] + (5.4) 

II. A^O,B> 4(^2 - K) >0, K <0: 

Eq. ( |4.9|) {A 7^ 0) leads to solitary wave ("bump" or "well" type solutions) for the 
function U [r]) in the following cases: 

2. U4 = U3 = U2 <U <Ui <0 

Eq. ( |5.5| ) represents an "algebraic bump" . 

3. U4<U <Us = U2 = Ui <0 

This is an "algebraic well" . 

4. U^<U3 = U2<U <Ui<0 

u = U2+ (u^-u^m-u,) 



{Ui - U^) cosh" A J {Ui - U2){U2 - U^){v - Vo) - {Ui - U2) 

(5.7) 



Eq. (|5.7|) is an "exponential bump" . 
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5. Ui<U <U3 = U2<Ui<0 



U = U' 



2 



(f/i - f/4) cosh^ AJ{Ui - U2){U2 - U,){r] - r/o) - (f/2 - U2 



(5.8) 



This is an "exponential well". 
Further elementary solutions of eq. ( [4.9|) are trigonometrically periodic. 
6. U4 = U3 < U2 < U < Ui < 0, K < 0, B > 4(^2 -K)>0 

U2-U3+ (f/i - U2) cos' AJ(Ui-U3)(U2-U:,)(,^ - ,,„) 



This type of solution also occurs only for the algebras dST 
7. f/4 < f/ < f/3 < f/2 = f/i 

(f/l-f/4)(f/l-f/3) 



U = Ui 

f/i - f/3 + (t/3 - U,) cos2 A J(f/i - f/4)(f/i - f/3)(r7 - r/o) 

(5.10) 

This solution can occur in the case of algebras ( |5.2|) for Ui < 0, i.e. all roots negative 
(and then conditions (5.5) hold). It can also occur for f/3 < < f/2 = f/i and this 
allows us to have i? < 0, Thus, solutions (|5.10|) can occur for all algebras (and variables 
^) 1-23 in Table |. Notice however, that they are periodic, rather than localized, in 
the variable rj. 

The remaining solutions are periodic and expressed in terms of Jacobi elliptic functions. 
We have: 

8. A = 0, > 0, B < -SK < 0, f/2 < f/ < f/i < 

U= e=^^l^ (511) 

U2 + (f/i - f/2) sn2 (v/-^(^-^o),A:) ' (-^^2) 



This occurs for the algebras (|5.3| ) 
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9. A = 0, K <0, B > -8K > 0, f/2 < t/i < f/ < 



i? = V^sny^(r^-r/o,fc), = (5.12) 



The algebras concerned are those of eq. ( |5.2| ) 

10. A ^ 0, f/4 < f/ < f/3 < f/2 < Ui 



{Us - f/4) sn^[l3ir] - r/o), k] + U,-Us ' 



(5.13) 



This can occur for Ui < 0, then we must have B > 4(yl^ + {~K)) > 0, i.e. the 
algebras ( |5.2|) . It can also occur for f/3 < f/4 < < f/2 < Ui, then all of the algebras 1- 
23 of Table H can occur. 

11. A^0,Ui<U3<U2<U<Ui<0 

_ U.jUi - U2) sn'lPiv - r^o), k] + [/^(f/^ - f/4) . 
(f/i-t/2)sn2[/3(r^-r^o),A;] + f/2-f/4 ' ^ ' ^ 



with k^ and /3 as in eq. ( ^.13] ). we must have B > 4(yl^ + {—K)) > and hence 
algebras 



12. A ^ 0, f/4 < f/ < f/i, f/2,3 =p±iq,q>0 

[CUi - DUi] cn[/3(r/ - r^o), A;] + Df/i + CU^ 



U 



iC-D)cn[P{ri-rio),k]+C + D 
C={Ui-pf + q\ D = {UA-pf + q^ (5-15) 



This situation can occur for all algebras 1-23 of Table III 



13. A = 0, K = 0, B < 0. We obtain the solution ( |4.17| ) for algebras ( p. 3D with as given 
in Table | (and r] given by eq. (t4.6|)). 
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The algebras No. 24-29 of Table || correspond to variables ^ such that (V^)^ = and 
hence to first order ODEs. The solutions are readily obtained and we just list them. 

No. 24: R = R^^ip = ay-bx+ ^ (x + t) + ^o, & ^ (5.16) 



No. 25: 



R = cV^t ± ^c\x + t)-l,^P = ^{t-x)- ^^^^ ^^1^ + ^0 (5.17) 



It^ — x"^ — 



No. 26: R = cVxTt±Jc^{x + t)-l,'^ = b\nJ — + ^Po (5.18) 

* \ X + 1 

b 

No. 27: R = Ro,tp = -ln{x + t) +ilJo (5.19) 

No. 28 and 29 provide nonconstant solutions only for 6 = 0. Then F{C,) and R{C,) are 
arbitrary functions of = x + t. 

All solutions presented so far are group invariant solutions in the standard sense of the 



words |p|-[TO|] 



Let us mention that the PDEs (|1.5| ) and (1.6) can be reduced to ODEs of the form (3.6) 
and (3.7), by the transformation ( p.l|) where ^ and a are any functions satisfying eq. ( p.3|) . 
The restriction is that p h, s and I must be functions of C,- Group theory generates solutions 
of these equations by the requirement that F and ^ in (p.l|) be invariants of subgroups of the 
symmetry group. However, other solutions may exist, corresponding e.g. to so called "null 
variables" PT| , P2[] , to "conditional symmetries" PB| , P^ , or simply generated by the "direct 
method" of Clarkson and Kruskal ||25|| . 

Let us just give some examples of such variables. 

First a few words about null variables and the corresponding solutions. Consider a 
variable ^ satisfying 

(VO^ = = (5.20) 



The equations for F{^) and R{^) reduce to eq. (3.13). As mentioned in Section |lTT| , if we 
have also 

(Va, VO = (Va)^ = Aa = (5.21) 
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(e.g. for a = const), then F{^) and R{C) are arbitrary functions. We have aheady encoun- 
tered this situation for C, = x + t, however eq. ( |5.20| ) have more general solutions pl| , |22| , 
that can be written in terms of Riemann invariants. 
Indeed let us put 

^ = H{a), a = (a, x) = a^t — aix — a2y, (a, a) = (5.22) 

where a is a lightlike vector, depending on ^ (i.e. eq. ( ^.22| ) defines ^ implicitly). It is easy to 
check that C, of eq. ( ^.22| ) satisfies eq. (eq5.20) for any choice of the function H and lightlike 



vector a(^). The function H{a) can be chosen to be H{cx) = a with no loss of generality, 
since and Ri^) are themselves arbitrary. Thus we can replace eq. ( p.22| ) by 



i = aoiOt - a^iOx - a2{i)y, a' = (5.23) 

Choosing a to be constant, we recover the variable ^ = x + t (up to a Lorentz transforma- 
tion). Other choices give different results, that become explicit if we can solve the algebraic 



equation ( |5.23| ). For example choose 



a=(l,e,Vl-a- (5.24) 



Solving eq. ( |5.23| ) for ^, we obtain 



'l + x)t±^{l+xy + y'^ -f^ 
[l + xf + y^ 



^ = ■ (5-25) 



Choosing F{^) and R{C) appropriately, e.g. F{!^) constant and R{C) with compact sup- 
port we obtain a localized solution (localized in the variable ^). 

An example of a "conditionally invariant" solution is obtained by putting 



Ix'^ + y"^ 



We have a{x, y,t) = and 



p = ^ = -l, g = - h = s = l = m = 0, B = (5.27) 

in eq. (3.6) and (3.7). These values could hence be added to those in Table |. 
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VI. COMMENTS AND CONCLUSIONS 



Inserting the variables ^ and a of Table | into the formulas of Section V we obtain a 
great variety of exact analytic solutions. 

Some of our solutions are (possibly upto phase factors, contained in the variable a) 
really solutions of the 1 + 1, or 2 + dimensional CP^ model. Thus, algebras 1, 11, 12, 16 
provide solutions depending essentially only on x and t. Similarly, algebras 2, 13, 14 provide 
essentially static solutions (independent of t). A sizable literature exists on static solutions 
[26-30]. Particularly interesting solutions of this type are obtained for algebra 2 when we 
have ^ = y/x'^ + y'^ and we take 6 = 0. We obtain elliptic function solutions ( |5.11|) and (5.16) 
as well as the elementary solution (4.17), or more explicitly 

W = iyo(\/a;2 + y2)V"'^ (6.1) 

where a = n is an integer (and is the azimuthal angle in the xy plane). This can be 
identified as an "n-soliton solution" (or instanton) and it has finite energy Our static 
solutions are not new: they are to be found e.g. among those obtained by Purkait and Ray, 
or earlier [26-30]. 

The same algebra gives rise to a very different type of solution. If we take b ^ 0, a = 
we express -R(0 terms of the Painleve transcendent Py, as in eq. (4.15). The phase is 
ip = bt + ipo s,o that we have a global rotation of spins in the horizontal plane. To our 
knowledge, this type of solution is new. 

Algebra 14 introduces a "helical" type variable ^. Solution (4.17) in this case is 

^ i?o(x^ + y2-jafe/2(l+a2)g-fe/l+a2arctany/Xgi6</./a_ j^g 2) 

This solution is multivalued, even for b/a integer. This type of variable and solution could 
be pertinent in condensed matter applications, concerning e.g. critical phenomena in mul- 
tilayered films. 

In general, our method provides us with "local solutions" , not necessarily defined for all 
of M^. The solutions are not necessarily single valued and they can have singularities for real 
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values of the variable ^. Moreover, in view of the existence of the light cone it is sometimes 
necessary to consider spacelike and timelike regions of space-time separately, since solutions 
in these regions may differ. Typical exemples of this phenomenon are provided by algebras 
17, 18 and 19. We list two variables ^ in Table | for each of these, one valid for t^—x^ —y'^ > 0, 
the other for x"^ + y"^ — > 0. In all cases we restrict to a = 0, in order to have m = in 



the table. The simplest solutions are given by eq. (4.17) for B < and { f)A\) for B > 0. 
In the case of algebra 17, the corresponding solutions are: 



W = tanh biVx ^' + y'-t- _ ^\ ^ ^2 _ ^2 _ ^2 < ^^'^^ 

2 \ X + t J 

The two solutions can be connected on the cone, however their derivatives will be dis- 
continuous in any case. 

Similarly, for algebra 18 of Table | we find the elementary solutions 



^/ — arctan ^ -x^ -y^ jy ^-ih arctanh xjt ^2 _ ^2 _ ^2 ^ g 

arctanh y/t ^2 ^ ^2 _ ^2 > 



W = tanh - 
2 



Vx^ + - 
arctanh 4o 

y 



(6.4) 



Finally, for algebra 19 of Table | we have 

W = tanh ^ (^arctan ^ ~ - ^o) e'^''^'^^^^'^^/^ x^ + y^ - ^ > 0, 

= y|/^g6/2arctanh-yt2^^x2^^/tg-i6arctany/x ^2 _ ^2 _ ^2 ^ g 

In many soliton-like problems in field theory we are interested in solutions which are 
regular in M^, i.e. which are valid at all times (though this condition is sometimes relaxed 
a bit) and which are defined for — oo < x < oo, — oo < y < oo. Among them particularly 
important are those whose energy is finite (as they describe localised "soliton-like" field 
structures). 

If we restrict our attention to such field configurations we see that we should consider 
the energy density for our fields. As the energy density is given by 

(6.6) 



\Wt 








[l + \W 


2J2 
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we see that this gives us 



P 



[1 + i?2]2 + [1 + ^2]2 



(6.7) 



where if) is given as in and R 



dR 



. We can rewrite {ifjf + ipl + ipl) = {a^ +al + ay) + 



F'^{Q + 'Cx + iy) + '^F{it(y.t + ixOix + ^ytty) and then substitute the expression for F given by 
(U) (when m = 0). 

To get the total energy we should integrate p over all space 



To perform this integration, in some cases, we can replace the integration over x and y by 
an integration over ^ and another conveniently chosen variable (which may have a finite or 



and an angle, while in the cases of 1, 4, 11, 12, 15 and 16 ^ involves only x and as our 
variables of integration we can use ^ and y. Clearly in these latter cases the total energy of 
any nontrivial solution is infinite. 

The most extreme case corresponds to the algebra 10. In this case ^ = t, energy density 



by 0'^ = (1 — -R^(^))/(l + -R^(O) aiid is independent of x and y, while 0^ and 0^ depend 
on X and y only through a. Thus treating 0* as components of a spin vector field (of unit 
length) we see that this solution describes very coherent movements of spins which move 
up and down in phase for all x and y and whose movements in the horizontal plane are 
modulated by a and F{t). 

Similar spin wave interpretations can be given to other solutions. In particular this is the 
case when the symmetry variable is more complicated than in the cases mentioned above. 
One can think of applications in condensed matter physics, the theory of nematic liquid 
crystals etc. and even in cosmology. In some of such systems the orientation of does not 
matter; such cases can be described by a larger class of our solutions. At the same time 
we can consider W{x,y,t) as a Landau-Ginzburg field which arises in many applications in 




(6.8) 



an infinite range). Thus in the cases of algebras 2, 3 



and 19 of the Table | we can use ^ 



is independent of x and y and so the total energy is infinite. In this case 0'^ of ( |1.3D is given 
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condensed matter physics (as can be checked the Landau-Ginzburg equation is very similar 
to the equation derived from ( [L.4| )). Indeed, at least one version of the Landau Ginzburg 
equation has been treated using the group theoretical techniques applied in this article. The 
context was that of magnetic phenomena in external fields. pT| 

However, returning to the field theory soliton-like applications, in which case the reduc- 
tions 2, 3, 17, 18 and 19 are particularly relevant, we note that using an angular variable 
of integration makes it more likely that a given solution will describe a time evolution of a 
field configuration of finite energy. 

Clearly it would be desirable to analyse further the physical implications of this and 
other solutions. We hope to be able to report on this in the near furture. 
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TABLES 



TABLE I. 



No 


Algebra 




a 


P 


9 


h 


S 


I 


m 


B 


1 


Ki + aZ, P2 + 6Z 




by — a arctanh j 


1 


1 

£ 













2(62^2 + „2) 


2 


Po + aZ, L + bZ 




bt + a arctan — 

X 


~? 


1 













2(62^2 _ a2) 


3 


Ki + aZ, K2 + L 


■\J t'^ — x^ — y'^ 


—a In \x + t\ 


-2 


1 

















4 


K2 + aZ, Po - P2 


X 


— a In \y + t\ 





1 

















5 


D + aZ, Po - Pi 


t+x 


aln |t + x\ 


-2 
{ 


1 

F 

















6 


Ki + eP2 + aZ, 
Po-Pi 


j/e + arcthf + iln(t2 _ 


aey 





1 








a2 








7 


D + e{K2 + L) + aZ, 

Po-Pi 




-aln|*±2| 





1 

















8 


D + aKi + bZ, 


2±*y(-i),a^l, 


hbHy) 


2-a 1 


^(2-a)/(a-l) 


6 1 


-b 1 


6^ 1 










a-1 J 


2(a-l) S 


2(a-l)^ 


4(a-l)^ e 




Po-Pi 




















9 


D + K^ + t{Po + Pi) + aZ, 
Po - Pi 


^x+ty-2c^ e = itl 


fix + t) 


l+2e 1 
2^€ 



















10 


Pi + aZ; P2 + bZ 


t 


by + ax 





1 








-(62+a2) 





2(62 +a2) 


11 


D + aZ, P2 


X 

t 


alnt 


2« 


1 


aS 
?^ 


a 







2a2 


12 


D + Ki+ e(Po + Pi) + aZ, Pa 


t + X — elii\t — x\ 


f ln|t-x| 





1 


ea 
4 











8 


13 


D + aZ, Pq 


arctan — 

X 


aln ■^x'^ + y'^ 





1 








a^ 





-2a2 


14 


D + aL + bZ, Po 


i ln{x^ +y'^) - ^ arctan ^ 


— arctan 

a X 





1 


i) 





6= 





,2 2 
r> a 


1+a^ 






15 


P2 + aZ, Po+bZ 


X 


by + at 





1 













2(a2 - 62) 


16 


D + aKi +bZ, P2, 


(t+xr+^ 


T^Jn\t-x\ 


1 


1 


b 













(t-x)i-» 

a > 0, ^ 1, 






25(a^-l) 


2(a^-l)^ 


17 


D + aZ, K2 + L + bZ 


(t2_^2_j^2)l/2 
x+t 


-b^+a\n\t + x\ 





1 


—a 





62 


a 








(x^+y^-t^)^/^ 


-bjf^+a\n\t + x\ 





1 


—a 

e. 





-62 


a 


2(1^+52) 


18 


D + aZ, Ki + bZ 


arctan (*^-^-^)^'^ 


-fearctanhf + f [t^ -x^l 





1 


_ g 
tan ^ 





62 -a2 


a 





arctanh 



(^2_^^2_,2)l/2 



-fcarctanhf + q\x^-t^\ 



tanh ^ 



sin 
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20 



D + aZ, L + 6Z 



D - Ki + aZ, 

Po - Pi + 6Z 

21 D + a/sTi + bZ, 

K2 + L 

22 D - + e(Po - Pi) + aZ, 

K2 + L 

23 D+ijCi+oZ, 
K2 + L + e(Po + Pi) 



arctan 



arctanh 



t 

(t2_^2_^2)l/2 



a ^ ±1, 
e = ±l 

6(t-x) + {t + x)^+6ty{t+x) 
[(t+a;)2+4ey]3/2 



-fcarctan | + § In -^a:^ + y2 q 
-6arctan| + §lny^x2+^ 
at- f (a; + t) + |ln|i; + t| -| 



tan g 

- g 
tanh ^ 

— 

2F 



■ln|a; + t| 



a\n\{t + xf + iey\ 



5 



1 5/6 



2 og 



1-a "''^•^ ^ 2 a+1 J ? a+1 2e " 



, i_£2 



■&2 



sin^ g 
g 

sinh"^ ^ 
g 

2F 

6(l-a) 1 

4(a+l)2 2(a+l)2 



^(l_a)/(l+a) 



4 g"' 2a 



8|l(l_^2)-l/3 



to 
00 



TABLE II. 



No Algebra ^ a (V^^ Aa (Va,VO (Va)^ 

24 P2 + aZ, Po - -Pi + 2; + 1 -6x + ay 6 -iP' - a? 

25 K2 + L + aZ,Po-Pi + bZ x + t bt- 6 6^ - |^ 

26 D + K^ + bZ,K, + L x + t iH^-x'-y') ..^l^^ 

27 - i^i + 6Z, + TT^^ |ln(t + x) -2^^^^+^ -b ^,.jj-/_\,y, 



28 D + iiTi + P2 x + t I ln(t - x) ^ 

29 D + iCi + Po - i'l + t a = 61ny 0^ 



29 



